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ABSTRACT 

We analyze the longitudinal 4-probe conductance of mesoscopic normal and su- 
perconducting wires and predict that in the superconducting case, large negative 
values can arise for both the weakly disordered and localized regimes. This con- 
trasts sharply with the behaviour of the longitudinal 4-probe conductance of normal 
wires, which in the localized limit is always exponentially small and positive. 
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In contrast with the huge literature on quantum transport in normal, phase 
coherent, structures |1], 0, a detailed knowledge of mesoscopic superconductoring 
structures is only now being developed. Such structures constitute new quantum 
objects, which during the past 3 years, have yielded many surprises || [|, [5]. Most 
theoretical work to date has focussed on two-probe transport coefficients, such as 
normal-superconducting (N-S) boundary || conductances and N-S-N 2-probe con- 
ductances |7|, ||. These are simpler to analyze than multi-probe coefficients and 
therefore without good reason for expecting new physics to emerge, there would 
seem to be little point in injecting unnecessary details into an already complex the- 
ory. The aim of this Letter is to demonstrate that multi-probe conductances contain 
new features which are absent from two-probe measurements. A key result of our 
work is that when a normal system with a positive 4-probe conductance is allowed 
to become superconducting, the sign of the conductance can change. 

As an example, we study the 4-probe conductances Gjk,i m of the 2 dimensional 
wire shown in figure la, where Gjki m = Ijk/(Vi ~ Kn) ; with VJ — V m the potential 
difference between voltage probes /, m and Ijk the current flowing from probe j to 
probe k. We consider the case for which the system size is smaller than the inelastic 
scattering length and focus attention on the longitudinal conductances (^13,24 and 
^14,23- If has been shown [[J that for normal materials, both negative and positive 
multiprobe conductances can both occur depending on the geometry of the device. 
On the other hand, numerical simulations of electron waveguide couplers [10], have 



found that for clean materials, the analogues of (^13,24 and Gi4,23 are always positive, 
in agreement with theoretical arguments valid for devices satisfying certain spatial 



symmetries 10, 11 . 



In this letter we first derive a general criterion which determines the sign of 
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the 4-probe conductance in any arrangement of normal conductors. In the case of 
localized normal wires, for which the quasi-particle transmission probability from 
one end of the sample to another is exponantially small, this shows that G13.24 
and (^14,23 are both positive and of order ~ exp — 2L/£, where £ is the localization 



length. In contrast, using a recently derived generalisation [12| of the multiprobe 
Biittiker formulae to the case in which Andreev scattering is permitted, we show 
that for a superconducting wire with exponentially small quasi-particle transmission, 
longitudinal conductances are finite and may be of either sign. We also present 
numerical results for the behaviour of 4-probe conductances in the presence of finite 
quasi-particle transmission along the wire. 

Consider first the case of a normal wire, for which it is known || that 

Gjk,im = D/ (TjiTkm — Tj m T k i) (1) 

where = J2 a /3 ^if an d T°f is the probability for a particle incident in channel (3 
of probe j to be transmitted to channel a of probe i. D is the determinant of the 
matrix obtained by crossing out any one row and column of the matrix of transport 
coefficients A, Aij = N5ij — Tij where N is the number of channels in each probe 
(assumed equal). The value of this determinant is independent of which row and 



column are removed |L3| , although because of the relation Ya Mj — Y,j Mj = the 
expression may be written in many equivalent ways. In order to make the analysis 
clearer we will write those scattering coefficients involving transmission along the 
wire, for example T 13 , T 24 (but not eg. T 12 or T 34 ) in the form Ty = exp — 2L/£ 
where has magnitude of order unity. We expand D by removing row 4 and column 
4 of A and substituting An = — Aij- This yields 

D — exp — 2-Zv/£ (7l2^23^~34 "I" ^~12^24^~34 "I" ^13^~21^~34 "I" ^14^21^34) 

+ exp -4L/f {T 12 t 2 ihi + T 12 t 2 4t 32 + T u t 13 t 23 + T u t 13 t 24 
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+721^14^31 + ^21^14^32 + ^34^14^23 + ^34^14^24) 

+ exp -6L/ £ (t 13 t 2A t 32 + t u t 23 t 31 + £14*24*31 + t 1A t 24 t 32 ) , (2) 

which demonstrates that for L » £, D decays exponentially with L. Although 
we have written D explicitly as a sum of powers of exp— 2L/£, expression (||) is 
exact and remains true when L <C £. The key observation here is that .D is the sum 
of positive terms and hence is always positive. Hence the sign of Gj^im depends 
only on the relative magnitudes of TjiT km and Tj m T k i. This was noted separately 



by Avishai and Band [i t] for a crossed wire arrangement and by Wang et al [10] for 
ballistically coupled wires. However both of these references rely on being able to 
apply certain symmetries to the system, which simplifies the form of D. The above 
analysis shows that D is positive, independent of such symmetries. 

Now consider the longitudinal conductances £13,24 and £14,23 in the limit L ^> £. 
Since none of T 12 , T 34 and T 43 are expected to decay with L, the denominator of 
expression (H) for these conductances is positive, with magnitude of order unity. 
Hence in this limit 

£13,24 ~ £14,23 ~ + exp -2L/Z (3) 
Since the two-probe conductance £2 of a device with substantial localization is 
of order ~ exp— 2L/£, we see that for normal wires in the localised limit 4-probe 
conductance measurements will give results of the same sign and order of magnitude 
as 2-probe measurements. 

Conductance formulae of the kind shown in equation (1), describe a normal 
scatterer connected to normal probes. The generalisation of this approach to the case 
where the scatterer incorporates superconductivity, but the probes remain normal, 
was first derived for 2 probes by 0, [8[ and more recently for many probes by 
Lambert, Hui and Robinson fl2|| . This generalisation leads to the introduction 
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of a matrix of transport coefficients A, which at zero temperature, has elements 
Aij = N5ij — T° + T-, where superscripts O and A refer respectively to normal and 
Andreev scattering. It is shown in that in general 

G ij)k i = dj (b ik - b jk - b u + b ki ) (4) 

Here, d = detv4, and b mn is the cofactor of the matrix element A mn . In order to 
highlight the behaviour of G in the localized limit we work to zero'th order in L/£ 
so that A becomes 
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(N-T°+T& -T°+T& \ 

-Tg+TA N-T& + T& 

N-T°+T£ -T°+T* 

V -Tg + Ti N-T& + T&J 



(5) 



The respective determinants dxi and dsR of the top left and bottom right blocks 
of A are both positive, because in general iV = Sj =1 (7|? + T-) for any i, so that 
An > A12, A 2 2 > A 2 i etc. Hence d = dTbdsR is positive in the localised limit. 

Now consider the longitudinal conductance £13,24 = dj (612 + 634 — &13 — ^24 ). To 
zeroth order in L/£ we have 613 = 6 2 4 = 0, 612 = —a2\dBR-, and 6 34 = —a i3 dTL- 
Hence 

1 21 1 21 _|_ 1 43 1 43 J ^ 
drL dsR ) 

A similar expression can be derived for £14,23 Since all the terms on the right hand 
side of OH) are of order unity, Gi3,24 would also be expected to be of this order. 
However, in contrast with the behaviour of normal systems the conductance can be 
of arbitrary sign, depending on the relative magnitudes of the normal and Andreev 
scattering coefficients for transmission from a current to a voltage probe at the same 
end of the wire. In general stronger Andreev scattering favours negative longitudinal 
conductances, while stronger normal scattering favours positive conductances. It 
should also be noted that it is possible in principle for the two terms in (|6|) to 
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approximately cancel, so that the longitudinal conductance may become arbitrarily 
large, even in the strongly localised limit [114 . 



To test whether finite negative multi-probe conductances can occur in practice 
we have performed numerical simulations of both normal and superconducting mul- 
tiprobe structures. Following Avishai and Pichard et al [|U|, |I6|| , we model a 2d wire 
with 3 channels and 4 probes by a network of Id wires connected at nodes, shown in 
figure lb. In the figure M denotes the number of slices between the pairs of probes. 
The figure shows the case M = 10, but simulations were performed over the range of 
lengths M = 1 to M = 40 for the normal wire and M = 1 to M = 100 for the super- 
conducting wire. Disorder is introduced by specifying random scattering matrices 
at each node, while allowing perfect transmission along each Id wire. For a normal 
system, the N w x N w scattering matrix of a node connecting N w wires is modelled by 
equating it to exp(iH), where H is a N w x N w Hermitian matrix chosen as follows: 
each element along or above the main diagonal is a real number chosen at random 
and independently of the other elements of H between — rr and tt, and the elements 
below the diagonal are chosen to ensure H = H T . For a superconducting sample 
the same procedure is adopted, except that the matrices are of size 2N W x 2N W and 
the elements of H are restricted to satisfy particle-hole symmetry. To obtain the 
scattering matrix for the whole network, we employ a numerical S'-matrix reduction 



algorithm [17, 18 1, details of which are explained more fully elsewhere [21 



Figure 2 shows the logarithm of the conductance as a function of the 

number of slices M in the wire. For each value of M, the conductances arising 
from 100 different realisations of the disorder are shown as dots in the figure. The 
inset shows corresponding results for the logarithm of the transmission coefficient 
T 13 . Clearly the typical values of both G1324 and T 13 decay exponentially with 
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M for large M. The different localization lengths for the two systems are to be 
expected, since the requirement that particle-hole symmetry must be satisfied for 
the superconductor will influence the level statistics of the random node scattering 
matrices. For a superconducting system, figure 3 shows corresponding results for 
^13,24 (plotted on a linear scale) and T13. This shows that for large M, whereas the 
transmission coefficient decays exponentially to zero, the conductance remains finite 
and can have arbitrary sign. This confirms our expectation in the strongly localized 
limit, based on equation (6) and additionally shows that negative conductances arise 
in the presence of quasi-particle transmission. 

In the above simulations positive and negative conductances occur with roughly 
equal frequency, because our choice of random node scattering matrices favours 
neither normal nor Andreev scattering. Historically, the phenomenon of Andreev 
scattering, which yields charge transport in the absence of quasi-particle transmis- 
sion, was used to explain the marked difference between thermal and electrical con- 
ductance across normal-superconducting interfaces. In the absence of quasi-particle 
transmission, the ends of the sample scatter quasi-particles independently and there- 
fore apart from a dependence of the condensate potential on the applied potential 
difference J7[ [EJ, the voltage probes become completely decoupled. In this limit, 
the voltage applied to probe 1 (3) serves to cancel the current due to quasi-particles 
from lead 2 (4) and therefore a 4-probe measurement can be viewed as two inde- 
pendent measurements of quasi-particle charge imbalance at the two ends of the 
system. In the absence of inelastic scattering there is no difference in principle be- 



tween charge imbalance measurements [22] and point contact spectroscopy of the 



kind described by Tsoi and Yakovlev p3"| . Here the sign of the voltage due to 
quasi-particles was shown to be reversed by the application of a magnetic field. By 
analogy such experiments, one expects negative 4-probe conductances of the kind 
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predicted in this letter, to be particularly sensitive to the presence of applied or 
internal magnetic fields. 

This work has benefitted from useful correspondance with Marcus Biittiker, who 
drew our attention to the results of ref. [13|]. Financial support in the UK from the 



Science and Engineering Research Council and the Ministry of Defence is gratefully 
acknowledged, and one of us (MJ) would like to thank the Science Foundations of 
the United States and Japan for a fellowship for research in Japan. 
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Figure Captions 

1. (a) A mesoscopic structure connected to four external probes numbered 1,2,3,4. 
l.(b) A network of one dimensional normal wires connected at nodes, with 3 wires 
in each external lead. The scatterer consists of M slices of nodes. For a normal 
system, the scattering matrix of each node scatters particles into particles. For 
a superconducting system, the scattering matrix of each node also incorporates 
Andreev scattering. 

2. The main figure is a scatter plot of the conductances Gi3,24 of a normal system 
as a function of the number of slices M. The inset is a corresponding plot of the 
transmission coefficients T 13 . All quantities are plotted on a logarithmic scale. 

3. Results for the conductances and transmission coefficients of a superconduct- 
ing system. In this case since the conductance no longer decays exponentially with 
M, only the transmission coefficients are plotted on a logarithmic scale. 



